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Abstract. The existence of true scale-invariance in slowly driven models of self- 
organized criticality without a conservation law, as forest-fires or earthquake automata, 
is scrutinized in this paper. By using three different levels of description - (i) a simple 
mean field, (ii) a more detailed mean-field description in terms of a (self-organized) 
branching processes, and (iii) a full stochastic representation in terms of a Langevin 
equation-, it is shown on general grounds that non-conserving dynamics does not lead 
to bona fide criticality. Contrarily to conserving systems, a parameter, which we term 
"re-charging" rate (e.g. the tree-growth rate in forest-fire models), needs to be fine- 
tuned in non-conserving systems to obtain criticality. In the infinite size limit, such a 
fine-tuning of the loading rate is easy to achieve, as it emerges by imposing a second 
separation of time-scales but, for any finite size, a precise tuning is required to achieve 
criticality and a coherent finite-size scaling picture. Using the approaches above, we 
shed light on the common mechanisms by which "apparent criticality" is observed in 
non-conserving systems, and explain in detail (both qualitatively and quantitatively) 
the diff'erence with respect to true criticality obtained in conserving systems. We 
propose to call this self-organized quasi- criticality (SOqC). Some of the reported results 
are already known and some of them are new. We hope the unified framework presented 
here helps to elucidate the confusing and contradictory literature in this field. In a 
second accompanying paper, we shall discuss the implications of the general results 
obtained here for models of neural avalanches in Neuroscience for which self-organized 
scale-invariance in the absence of conservation has been claimed. 
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1. Introduction: Critical Self-organization with and without conservation 

Power-law distributions are quite common in Nature: earthquakes and micro-fractures, 
solar flares, weather records, snow avalanches, crackling, and 1/f noise are just a 
few examples of systems displaying scale-invariance [1]. Many of these, among many 
others, have been claimed to be critical, i.e. to lie at, or very close to, a critical 
point, which ensues scale-invariance and the concomitant power-law distributions. It is 
worth stressing that power-laws (or approximate power-laws) emerging from the complex 
interactions of many- component systems placed at the vicinity of a critical point and, 
therefore, with diverging correlation lengths and the associated power-law decay of 
temporal and spatial correlations, should be clearly distinguished from other power-law 
distributions arising in many different contexts (word distributions, city populations, 
citations, etc). These latter can be generated by a wealth of non- critical mechanisms 
as, for instance, multiplicative noise or fragmentation processes (see [21 [3l H] for recent 
reviews), without the need to invoke criticality. 

For a system to be critical, it does not suffice to have power-law distributed 
observables but, more crucially, it has to obey finite-size scaling: measurements 
at various system-sizes (scales) can be related to each other by re-scaling variables 
and quantities in some specific "scale-invariant" way (i.e. scaling collapses can be 
performed) . 

Given that, in standard phase transitions (both in equilibrium and away from 
it), a precise parameter tuning is required to reach criticality and generate power- law 
distributed quantities, an alternative explanation for the emergence of generic critical 
scale-invariance (i.e. criticality occurring without requiring fine-tuning) was historically 
much needed [5]. How does criticality emerge spontaneously? 

In a seminal paper, Bak, Tang and Wiesenfeld (BTW) [6] introduced, back in 1987, 
the concept of self-organized criticality (SOC) [3, El El [IDl IHl US] aimed at solving the 
previous conundrum. The research line opened by their breakthrough work continues 
to attract, more than twenty years after, a great deal of interest. From this perspective, 
the overwhelming activity generated in the last decade around scale-free networks [13] 
can be certainly considered as a prominent extension of preceding research on SOC. 

A handful of mechanisms were proposed under the common name of SOC to justify 
the abundant presence of scale-invariance in the natural world. The most successful 
among them is the one exemplified by the original BTW sandpile model. Many 
variations of the BTW sandpile were proposed: the Manna sandpile [H], the Oslo 
rice-pile model [15j, and the Zhang model [16], are just a few of them. For the 
sake of completeness and for future reference, a brief description of the best-known 
prototypical SOC models is provided in Appendix A. Other possible mechanisms, as 
extremal dynamics [H], have been studied, but we shall not be concerned with them 
here. 

The idea inspiring sandpiles, ricepiles, and related SOC models is that many 
systems in Nature, when pushed/driven slowly, respond very irregularly with rapid 
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re-arrangements (avalanches, bursts) of a broad variety of sizes. The distribution of 
such avalanches is scale-invariant in many cases. Of course, sandpile SOC models 
are too simplistic to reproduce the detailed behavior of real sand-piles; they can be 
regarded as "metaphors" capturing only some particular features of real systems in 
a stylized manner. Nevertheless, after many (partially) failed attempts [H] power- 
laws for avalanche in real granular piles, obeying finite-size scaling, were experimentally 
measured [19]. Moreover, other physical situations, as vortex avalanches in type-II 
superconductors, can be mimicked as sandpiles, and their critical properties rationalized 
in terms of these [20] . 

The dynamics of a generic sandpile model can be synthesized as follows: some 
type of "energy" or "stress" (sand-grains) is progressively injected in discrete units 
(i.e. grains are dropped singly) at the sites of a spatially extended system (usually a 
two-dimensional square lattice) at a slow timescale. Whenever a certain threshold of 
local energy is overcome, the corresponding site becomes unstable and its accumulated 
energy is redistributed at a much faster timescale among its neighbor sites. These, 
on their turn, can become unstable, and trigger a cascade of re-arrangements, i.e. an 
avalanche or outburst of activity. Local redistribution rules are conserving in sandpiles: 
energy does not dissappear, but only diffuses around. Open boundaries are customarily 
considered to allow for energy release from the system. Once all activity ceases (i.e. 
the avalanche stops) new energy is injected (i.e. a grain is dropped) into the system, 
and so on, until a statistically stationary state is reached. In such a steady state 
avalanches are scale invariant, i.e their sizes/times are power-law distributed up to a 
maximum scale imposed by the system size: the system "self-organizes" to a critical 
point [3 [HI [9l [TOl [m [2H [22l [23] . Moreover, the associated power spectrum exhibit 1//° 
noise [21]. 

For the sake of generality, let us fix a nomenclature common to all models discussed 
in this paper: "energy" refers to the accumulated and transported magnitude and 
"activity" describes "energy above threshold". In some cases, sites below threshold 
will be sub-classified in two groups: critical, which can become active upon receiving 
one input of energy (for instance, one grain) and stable, which cannot [25] . 

An essential ingredient of SOC is slow driving (driving and dynamics operating 
at two infinitely separated timescales [TJ [HI [SI [211 EE] , i-e- avalanches are instantaneous 
relative to the timescale of driving). Such an infinite separation is usually achieved 
by driving the system only when all activity has stopped, but not during avalanches. 
For any finite separation of timescales, a finite characteristic (time/size) scale appears; 
hence, slow driving is a crucial requirement for generic scale-invariance to emerge 
[211 [23 [281 [29]. 

It was soon emphasized that energy conservation is also a key player for criticality 
to emerge in sandpile models [2ll [301 IMJ- Note that by "conserving system" one can 
refer either to models with bulk conserving dynamics and boundary dissipation, or to 
cases with a bulk dissipation rate vanishing in the large system size limit [32]. Various 
relatively simple arguments were proposed to rationalize the existence of true criticality 
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in the steady state of conserving self-organized systems (some of them are briefly 
discussed in the next section). However, these arguments cannot be easily extended 
to similar non- conserving systems. In particular, different works have shown that the 
level of dissipation acts as a relevant parameter in the renormalization group sense: any 
degree of bulk dissipation breaks criticality in sandpile models [211 [3ll 133] . 

Still, given the large variety of natural phenomena exhibiting (exact or approximate) 
scale-invariance in which some form of dissipation is inevitably present (i.e. systems 
without any obvious conserved quantity), alternative mechanisms for self-organization 
to criticality in the absence of conservation were needed to achieve a comprehensive 
picture of generic scale-invariance pi] . 

Two acclaimed non-conserving self-organized models, or better, two families of 
models proposed to fill the gap between theoretical understanding and empirical facts 
are earthquake and forest-fire models (see Appendix A for definitions). These are highly 
non-trivial, and interesting models with rich and complex phenomenology. Owing to the 
lack of solid theoretical arguments, analogous to the ones sketched above for conserving 
systems, and despite of numerical evidence showing power-laws for some decades, the 
existence of true generic scale-invariance in them has been long controversial. 

It is beyond the scope of this paper to review exhaustively the large body of 
interesting literature devoted to non-conserving SOC models, some aspects of which 
remain unsettled. But, let us just underline that the state-of-the-art is, as documented 
in the next section, that none of the considered non- conserving models is truly critical; 
they just show "apparent scale-invariance" or "dirty criticality" for some decades. 

However, this final conclusion has not been sufficiently stressed and it has certainly 
not permeated the literature. This is likely due to the absence of a general theory, 
which may suggest that the results discussed above are specific to each particular 
model. Indeed, works continue to be published assuming or claiming true criticality 
for SOC non- conserving systems. For instance, in [M] an interesting and solvable "non- 
conserving model of SOC" was proposed and studied analytically. In a more recent series 
of papers, Juanico and collaborators claim to have constructed different non-conserving 
self-organizing models with applications in various fields (as neuroscience, population 
dynamics, etc.) |35j- Also, in a recent work, Levina et al. propose a non- conserving SOC 
model for neural avalanches to capture the apparent scale-free behavior of avalanches 
of activity observed experimentally in networks of cortical neurons |36]. An exhaustive 
analysis of this last model, as well as a study of the possible relation between SOC and 
neural avalanches, is left for a separate publication. 

The aim of the present paper is to put together some previously existing results, 
scattered in the literature, (although this is not intended to be an exhaustive review 
article) and, more importantly, to rationalize the conclusion that none of the above 
mentioned non-conserving models, nor variations of them, exhibits true criticality, 
within a unified framework. To this end, we rely on different types of analytical 
arguments complemented by computer simulations. In passing, we shall report on a 
number of new results and present a critical discussion on the existence of true scale- 
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invariance in Nature. 

The rest of the paper is structured as follows: in Section [2|, we briefly review 
conserving and non- conserving models of SOC, as well as some arguments to justify 
the existence of criticality in the first group. In the remaining Sections, we elucidate 
the existence or not of criticality in non-conserving systems using different approaches 
of increasing complexity. In particular, in Section [3l we discuss a simple mean-field 
approach based on an energy balance equation; it is useful to illustrate some key 
concepts as the "loading mechanism" . In Section HI we study a self-consistent mean- 
field approximation, namely the so-called self- organized branching process; it serves as 
an adequate benchmark to scrutinize the effects of dissipation and "loading" in critical 
self-organization. In Section [5], we present (and briefly review) the Langevin theory 
of conserving SOC systems. It constitutes a solid basis to implement dissipation and 
loading in a systematic way and to provide clear evidence on the lack of criticality in non- 
conserving models. Finally, the conclusions and a critical discussion of the implications 
of our main results are presented in Section El 

2. Conserving versus non-conserving models of SOC 
2.1. Conservation and criticality. 

As said above, different type of arguments of different nature justify the existence of 
true criticality in conserving SOC models. Some of them are as follows: 

• The energy, introduced into a pile at generic sites, can reach the boundaries 
(and, thus, be dissipated) only by means of the diffusive transport of grains 
occurring during avalanches. Owing to this, and provided that a steady state exists, 
arbitrarily large avalanches (of all possible sizes) should exist for an arbitrarily 
large system size, ensuing a power-law size-distribution. Contrarily, in the presence 
of non-vanishing bulk-dissipation, energy disappears at some finite rate, and 
avalanches stop after some characteristic lifetime/size determined by the dissipation 
rate l2ij. 

This type of argument, even if commonly used in the literature, is (at best) 
incomplete, and can be misleading. It does not consider the possibility of having 
a characteristic scale larger than the system size, which would allow for avalanches 
to reach the boundaries. Actually, this is what happens in many real sandpiles 
(with inertial effects): energy is dissipated quasi-periodically in large system- 
wide avalanches rather than in scale-invariant avalanches (see the chapter on 
experimental set-ups of sandpiles in [HI El [37] or [21]). 

• From a more abstract viewpoint (not referring specifically to sandpiles or SOC), 
energy conservation follows from the existence of a continuous symmetry (in this 
case, temporal translational invariance) as a consequence of Noether's theorem [38]|. 
This also holds the other way around: time translational invariance implies energy 
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conservation. When energy conservation is violated, the corresponding symmetry 
is broken and a characteristic (finite) time-scale appears generically. 

• From a field-theoretical perspective, in order to have scale-invariance, generic 
infrared divergences are required. But these are generically lost in the presence 
of a non- vanishing linear "mass" term (adopting the field-theory jargon) as, for 
instance, a dissipative term. In particular, if a term —eip, is introduced into the 
simplest mesoscopic equation for a diffusive field 'ip{x,t) 



where ri{x, t) is a zero-mean Gaussian white noise, a simple calculation reveals that 
the equal time two-point correlation function can be written as 



for X much larger than the correlation length ^ ~ 1/ i/e. Accordingly, it is only 
for e = that ^ diverges, the exponential cut-off disappears, and the correlation 
function decays algebraically; for any non-vanishing value of e there is a size- 
independent exponential cut-off. Something similar occurs for other correlation 
functions. 

Note that the noise in Eq.([T]) is not fully conserving, but only conserving on average 
(i.e. conservation needs only to hold on average to preserve scale invariance |39j). 
The same conclusions can be also deduced for an equation analogous to Eq. ([1]) but 
with a strictly conserving noise (see [211 123, 128]). 
• Last but not least, a mesoscopic Langevin equation that captures the critical 
properties of stochastic sandpiles and related models with a conservation law has 
been proposed [221 SOI El 112]. It describes systems with many absorbing states 
(which correspond to the many stable microscopic configurations of a sandpile) and 
a conservation law. While a detailed description of this is left for a forthcoming 
section, we just stress here that it constitutes a sound field theoretical representation 
of conserving SOC, reproducing all critical exponents. The underlying mechanism 
of SOC, highlighted by this theory, cannot be straightforwardly generalized to non- 
conserving systems without including an additional fine tuning (as we shall show 
in Section [5|). 

In summary, there exists solid theoretical ground to underpin the existence of true 
criticality in conserving self- organized systems. The same type of arguments cannot 
be easily extended to non-conserving systems. In particular, as said already, non- 
conservative sandpiles have been explicitly shown to be non-critical. 

2.2. Critical self-organization without conservation? 

The two main prototypical non-conserving "self-organized" models (or families of 
models), studied profusely in the literature, are (see Appendix A): 



dtij{x, t) = D^V^ij{x, t) + r/(f , t) 



(1) 



(V^(x,t)v^(o,t)) 



exp(-x/0 



(2) 
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• Earthquake models as the Olami-Feder-Christensen (OFC) [S] and related ones 
[HI HHl im SHI So], in which the degree of dissipation can be explicitely controlled, 
and for which conservation is preserved only for some specific choice of a parameter. 

• Forest fire models as the Drossel-Schwabl automaton [50] (and earlier versions of 
it [51j), in which there is no conserved quantity whatsoever. 

Together with bulk-dissipation, these two models have a common key ingredient, 
absent in conserving systems: there is an increase of the "background energy" at some 
(or at all) sites, occurring between avalanches; (i) the accumulated stress at each site 
grows continuously between quakes in earthquake models as the OFC, and (ii) new trees 
grow between two consecutive fires in forest-fire models (see below). 

The effect of these "loading mechanisms" , as we shall call them generically, is to 
counterbalance the loss of "energy" (grains, stress, trees) produced by dissipation and, 
in this manner, try to restore conservation on average and, thus, criticality [25l |3H |35] . 
However, let us caution that such a compensation needs to be exact and, therefore, 
unless a new mechanism giving rise to a perfect cancellation is devised, fine-tuning of 
the loading rate is the only obvious way in which conservation can be restored. 

Let us now discuss in more detail these non-conserving archetypical models. 

2.2.1. Earthquakes: The Olami-Feder-Christensen cellular automaton [44] is a 
simplified version of a previously proposed fault dynamics models: the spring-block of 
Burridge-Knopoff model [15] and related stick-slip models [ITJ HH] designed to capture 
the essence of earthquakes (the Gutenberg-Richter law [52] for the distribution of 
magnitudes) as well as of similar systems with friction and jerky motion (see Appendix 
A). 

At each time step, the "forces" F{i,j) (or energies, to stick to our generic 
terminology), defined at each site of a two-dimensional lattice, are increased at a 
constant rate. Whenever F at any site reaches the threshold value, Fthr = 1, it is reset to 
zero and the forces at its 4 nearest neighbors are increased by an amount 6F = aF{i,j). 
This might trigger cascades of re- arrangements, i.e. avalanches. Observe that the bulk 
dynamics is conserving only for a = 1/4 in two- dimensions. 

Early computer simulations and theoretical results [HI [531 El] seemed to support 
the existence of criticality for values of a as low as a ~ 0.05. It was also early 
reported that, imposing periodic boundary conditions, the OFC model enters a cycle 
of periodic configurations with no sign of criticality whatsoever [55| l28l [56] . This 
suggests that the bulk dynamics is profoundly influenced by boundary conditions. 
Actually, it was proposed that boundary-induced heterogeneity is essential to obtain 
partial synchronization between clusters of different sizes and that such a partial 
synchronization or "phase- locking" mechanism is at the basis of the OFC complex 
behavior [55l EHl EE] . 

The role of different features (as, for instance, changes in the boundary conditions, 
introduction of quenched disorder in the local rules, lattice topology, etc.) on 
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synchronization and their effects in the properties of the OFC model have been largely 
analyzed in the literature [53 Ell Ell [53 |60l ill |62l |63l [M]. It has been also shown 
that results are affected by numerical precision [70]. The overall picture is that 
the "synchronization mechanism", even if fascinating, is too fragile as to be a solid 
explanation for generic emergence of criticality. A nice and rather exhaustive review of 
the literature on the OFC model and variations of it can be found in [72] . 

On the analytical side. Broker and Grassberger [65] and Chabanol and Hakim [66], 
in two independent papers, were able to calculate the energy distribution, the effective 
branching ratio, and the average avalanche size for a random-neighbor version of the 
OFC model, which turns out to be analytically solvable. Their main conclusion is that 
it is only in the conserving limit that the model becomes critical, while exponential 
cut-offs appear for any a ^ 1/4. Similarly, de Carvalho and Prado claimed, relying on 
an effective branching ratio analysis [67J, that the OFC model is only critical in the 
conserving limit (see also [68]). 

Remarkably, it is also shown in [65] and [66] that the average avalanche size is 
distributed as a power-law with a cut-off function, exp {jz^-, which diverges in a very 
fast way when the conserving limit is approached. This provides an explanation for the 
relative large power-law regimes observed even in the non-conserving case. It would be 
certainly nice to have extensions of this result to other non-mean-field like systems. 

In a similar line of reasoning, Kinouchi and Prado introduced the concept of ^^robust 
criticality" or almost criticality" [69]: for a fixed dissipation rate, systems with a loading 
mechanism are closer to criticality than systems without it. The reason for this is simple: 
moderate loading partially compensates energy dissipation. 

Finally, the most recent and exhaustive analyses by Miller and Boulter |71j . 
Grassberger [55], and Drossel and coauthors [72] conclude unambiguously, using a 
variety of arguments and large-scale computer simulations, that the spatially extended 
version of the non-conserving OFC model is not critical. 

As a consequence, the state-of-the-art is that, despite of the apparent power law 
distributions spanning for a few decades, the OFC model is not truly scale-invariant, 
except for its conservative limit. The question of whether real earthquakes are described 
or not by this type of SOC models or other type of mechanisms need to be invoked 
remains unsolved [731 EH [75] . 

2.2.2. Forest fires: The Drossel-Schwabl model [50] (see also [76]) is an improved 
version of an older forest-fire model [51] proposed to explain the apparent scale- 
invariance of real forest-fires [77] (see Appendix A). 

Three type of states are defined: z = or empty, z = 1 or occupied by a tree, 
and 2; = 2 or burning tree. At each time step, new trees grow, at rate p, at randomly 
chosen sites provided they were empty, and trees catch fire at a much smaller rate /. 
Fire propagates deterministically to neighboring occupied sites and, after burning, trees 
become empty sites. The relevant parameter is 1/6 = p/ f ^S\, and the model has been 
claimed to be critical provided that the double limit / — )■ 0, p — with //p — )■ 0) is 
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Observe that a double separation of time scales is imposed in the model definition: 
trees are born at a much faster rhythm than fires occur and fires propagate at a much 
faster pace than trees grow [50]. This is to be compared with the single time scale 
separation in sandpiles [50l El]. We shall discuss later the consequences of such a 
double separation of timescales. 

Analytical results and mappings into a branching process [79] first suggested some 
similarities with standard percolation models [SDl EOl EI]- Given the limited analytical 
tractability of these models, the controversy about the existence of true criticality was 
mainly played on the ground of computer simulations [50l l82[ l8Tl [83] . For sufficiently 
large systems, anomalies were reported to appear; among them: (i) the repulsive 
character of the fixed point {f /p = 0), (ii) the coexistence of largely-sub critical and 
supercritical clusters of trees, (Hi) the existence of two length scales with different 
exponents into the system, (iv) the violation of standard scaling for the distribution 
of avalanche sizes P{s), and (v) a pathological finite-size behavior [50, EH ESI ES, El]- 

Finally, when "massive" simulations of extremely large systems very close to the 
critical regime were accessible [851 EE], these anomalies turned into a lack of true critical 
behavior, beside of the apparent scaling observed for a few decades: the self-organized 
stationary state of the Dross el- Schwab I model is not critical. 

In the rest of the paper, we shall rationalize and generalize the above conclusion 
(i.e. absence of bona-fide criticality of earthquake and forest-fire SOC models) to generic 
non- conserving systems, to this end, we shall employ three different unified frameworks, 
as described in the three forthcoming sections. 

3. A simple mean-field approach 

As already mentioned, the random-neighbor version of the OFC model has been solved 
analytically, with the conclusion that, except for the conserving limit, it is not critical, 
but generically subcritical [651 ES] . On the other hand, in a subsequent work, Pruessner 
and Jensen [M] considered a modified version of such a model in which, by including 
a different (stronger) loading mechanism, they showed that criticality can be restored 
in the infinite size limit. In this section, we review the results in [31] and study the 
finite-size scaling of this and related mean-field systems. 

The model in [31] is somewhere in between forest-fire and earthquake models. It is 
defined as follows: consider a set of sites, each of them with an associated energy Zi 
(with z a continuous variable). As in the OFC model, three types of states exist: stable, 
with an energy < Zi < 1 — a; susceptible, with 1 — a < Zi < 1; and active sites, with 
Zi > 1. The main difference with respect to the OFC model is that, between avalanches, 
driving and loading operate as independent mechanisms [M] : 

• Background loading: in the spirit of forest-fire models [501 EH], (1/^) sites are 
randomly chosen and their respective energies are increased to 1 — a (i.e. become 
susceptible) provided they were stable; otherwise, nothing happens. 
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• Triggering of an avalanche: A randomly chosen site, i, is activated {zi = 1) provided 
it was susceptible. 

The relaxation dynamics within avalanches is identical to that of the random 
neighbor OFC model with m random neighbors: sites above threshold are emptied, 
z = 0, and the energy of its (randomly chosen) neighbors is increased by a fixed amount 
a. Conservation holds for a = 1/m. 

At a mean-field level, the condition for stationarity is given by the following energy- 
balance equation [M] : 



where the different terms are as follows: 

(i) l.h.s: For each relaxation event at site i, the amount of dissipated energy is 
(1 — ma)zi] hence, the average dissipation during an avalanche is (1 — ma)za{s), where 
(s) is the average avalanche size and Za the average energy of active sites. 

(a) r.h.s. first term: Triggering increases the energy of the selected susceptible site, 
i, by an amount 1 — Zi. The corresponding average increase is 1 — Zg, where Zg is the 
average energy of susceptible sites. 

(Hi) r.h.s., second term: Every time the background is loaded, the energy of a 
stable site, i, is increased by an amount [1 — a — Zi). This is attempted (1/6*) times 
and, the average number of triggering events before an avalanche is actually generated is 
(st/Cs (where (st and Cs are the density of stable and susceptible sites, respectively). The 
average increase of background energy per avalanche is, finally, {l/9){l — a — Zst){(st/Cs), 
where Zgt stands for the average energy of stable sites. 

In this way, Eq.(l3]) establishes that, for a steady state to exist, the average dissipated 
energy should be compensated by the averaged energy increase of driving and loading. 
Now, imposing in Eq.([3]) that (1/6*) diverges, one of the following two conditions must 
be obeyed for Eq.(^ to hold [3l]: 

(i) (1 — ma) = 0, i.e. there is strict conservation, or 

(ii) (s) diverges, which is a necessary condition for criticality. 

In the second case, by studying the probability distribution function of z and using a 
mapping into a branching-process, it has been shown that not only (s) diverges, but 
also that the system is critical in the infinite size limit [31]. Moreover, as expected for 
a mean-field model, the size- avalanche exponent is found tober = 3/2 [HIHIIHI 125] . 

However, as already pointed out in [M], for any finite system size, N, neither 1/6 
nor (s) are infinite. In such finite value of the parameter 1/6 must be fine-tuned 

to some precise value, l/6c{N), for Eq.(|3]) to hold. Such a value should diverge slower 
than N, in order to achieve the right limit {l/6c{N))/N — t- for — oo, but there is 
no analytical prescription in [3l] on how to fix it for each system size. 

This is a well-known problem, shared by forest-fire models, where the number of 
trees grown between two fires (1/^) is a parameter which needs to be carefully tuned for 





Self- organization without conservation : true or just apparent scale-invariance? 11 




10° 10^ S 10" 10° 10^ S 10* 



Figure 1. Left: Avalanche size distribution for the random- neighbor model 
introduced by Pruessner and Jensen [34], for = 2^^ and dissipation parameter 
a = 0.15. Subcritical, almost-critical (but still cut-off by system-size) and supercritical 
distributions, corresponding to different values of the loading parameter (1/^? = 2, 8, 
and 22, respectively) are shown. Right: Avalanche size distribution at the fine-tuned 
critical values of l/9c{N) for different sizes (from 2^^ to 2^^); the corresponding cut- 
offs grow with system size in a scale-invariant way. In the inset, the same data are 
collapsed into a unique scaling curve by using Eq.([S]). 



any finite size: too small values lead to subcritical fires, while too large values generate 
super-critical fires spanning the whole system (and generating a bump for large values 
in the size distribution). 

This is graphically illustrated in Fig. [1] (left), where the avalanche size distribution 
for the Pruessner- Jensen model is plotted for a system with = 2^^ sites, a = 0.15, 
and three different values of 1/6{N)] even if the values of 1/9{N) are large, the size- 
distributions are not pure power-laws: they are either subcritical (with an exponential 
cut-off) or supercritical (with a bump for large avalanches). Nevertheless, partial scaling 
is observed in any case. Given that the control parameter 1/6 is an integer number, 
criticality cannot be tuned with arbitrary precision (specially for small system sizes) but, 
still, for each value of A^ it is possible to find an almost critical value of 1/6, l/6c{N). 
In Fig. [1] (right) we show the avalanche size distribution at the fine-tuned critical point 
for different values of A^. Increasing the system size we have observed that such a value 
scales as l/6c{N) ~ A^0-33{5) suggesting 

l/6c{N) ^ N^/^. (4) 

A collapse of the critical size-distribution curves for different values of A^ (see the inset 
of Fig. H]) leads to P(s, A^) ~ s'^-^^^) exp(-s/A^°-6^(^)), compatible with 

P(s, N) ~ exp{-s/N^/^), (5) 

where the mean-field exponent r = 3/2 is recovered. 

In Section |5l we shall introduce a scaling Langevin theory for non- conserving SOC 
models, which explains in a straightforward way the two (formerly unknown) scaling 
laws, Eq.@ and Eq.l^. 
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Let us comment on the peculiarity of the thermodynamic hmit in this model (as 
well as in the Drossel-Schwabl forest-fire): the condition l/6c{N) — )■ oo is automatically 
fulfilled by imposing the double separation of time-scale discussed above, i.e. the 
second separation of timescales is tantamount to fine-tuning 6c{oo) to its critical value 
6c{oo) = 0. Note also that the infinite size limit is somehow pathological as the entire 
super-critical phase (as well as the critical point itself) collapses into a unique single 
point Oc{oo) = 0. Instead, for finite systems, a precise (not infinite) double separation of 
scales is required to set the system to the critical point, separating distinct sub-critical 
and super-critical phases. This boils down to the need of fine-tuning l/6c{N) for each 
value of to have a coherent finite-size scaling. 

On the contrary, in the conserving limit, large avalanches spanning the whole system 
are observed for any size and criticality is reached without resorting to careful tuning. 

Summing up, even if the random-neighbor model studied in [34| exhibits infinite 
avalanches and is critical in the infinite system size limit, it lacks of a well-defined 
finite size scaling and, therefore, it is not truly scale-invariant: for any finite system, 
deviations from criticality are observed if the control parameter 1/6{N) is not fine tuned 
to a precise iV-dependent critical value. In conclusion, this model does not qualify as a 
bona fide self-organized critical system. 

In this respect, the situation for the RN-OFC model studied in [651 [66] is even worse: 
given that it lacks a parameter analogous to 1/6' to be tuned, the degree of loading cannot 
be regulated and the model is generically subcritical for any non-vanishing dissipation 
rate even in the large system size limit [87] . 

4. Self-organized mean-field approach: Self-Organized Branching Process 

In this section, we complement the mean-field approach of the previous one by exploiting 
the self- organized branching process introduced by Zapperi, Lauritsen, and Stanley in 
[88] . First, in Subsection 14. H we introduce this approach for a conserving sandpile model 
(the Manna cellular automaton). Then, following also Zapperi et ai, in Subsect. 14.21 
we move on to analyze dissipative models, showing that they are generically subcritical. 
Finally, as a last step, in Subsect. 14.31 we implement a loading mechanism in the self- 
organized branching process which captures the essence of earthquake and forest-fire 
models, and explore under which circumstances the resulting model is critical. 

4.1. Conserving case 

In sufficiently high spatial dimensions (i.e. in the mean-field regime), avalanches in 
sandpiles rarely visit twice the same site; activity patterns are mostly tree-like. An 
avalanche can be seen as a branching process |79] in which an individual [ancestor) 
creates a fixed number k of descendants with probability p{k). The average number 
of descendants per ancestor, a = "^f^kplk), is called branching ratio. For a > 1 
avalanches propagate indefinitely (super-critical phase), for a < 1 they stop after a 
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Figure 2. Left: Sketch of the propagation of an avalanche in a self-organized branching 
process of size N = 2'^+^ — 1 = 15. The avalanche has size s = 7 (red spots), and has 
reached the border (i.e. there are two occupied sites in the last generation; <;{p, n) =2). 
Figure adapted from [88] 

typical number of generations (sub-critical phase), while the process is critical in the 
marginal case, Uc = 1 [79] . 

In this static branching process, fixing, without loss of generality, the number of 
descendants to /c = 2, a given active site branches in 2 with probability p or has no 
offspring with probability 1—p (see Fig. [2]), ensuing a = 2p and a critical value Pc = 1/2. 

Let us consider, as a simple example, the Manna sandpile with critical threshold 
Zthr = 2 (see Appendix A), and map its mean- field version into a self- organized branching 
process, in which p itself is a dynamical variable |88] . In the Manna dynamics, each grain 
arriving at site i can either generate activity (energy above threshold) if previously 
Zi = 1, or not if Zi = 0. An avalanche in high spatial dimensions can be, therefore, 
seen as a branching process with p = P{z = 1), i.e. the background energy density is 
nothing but the branching ratio [H9]. A generation is defined as the set of sites probed for 
activation at each time-step; after m generations there are N = 2™+^ — 1 involved sites 
(see Fig. 12]). For computer simulations, we fix a maximum number of generations m and 
impose (to mimic boundary dissipation) that at the m-th generation all grains are lost. 
Note that, apart from such a dissipative boundary, the bulk dynamics is conserving. 

Each avalanche modifies the background in which the next avalanche is to be 
started; i.e. it changes the value of P{z = 1). Thus, the branching probability p 
becomes a fiuctuating variable, as illustrated in Fig. [3l In the left part of the figure, the 
value of p is plot as a function of the avalanche number for different system sizes, while 
in the right figure the statistically stationary distribution of values of p is represented 
for various sizes; the width of the distributions decreases with increasing size and can 
be made as small as wanted. 

To recover analytically these computational observations, let Z{n) be the total 
number of grains into the system after n avalanches; then, p{n) = Z(n)/N. If q{p,n) is 
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Figure 3. Left: Evolution of tlic brancliing probability p{n) in a computer simulation 
of the self-organized branching process as a function of the number of avalanches per 
site, for different sizes N (ranging from 2^ — 1 to 2^^ — 1). In all cases, p fluctuates 
around p* = 1/2; the larger the system, the smaller the amplitude of fluctuations. 
Right; Stationary probability distribution functions of values of p for various sizes 
(same color code in both figures) ; the width decreases with increasing system-size and 
vanishes in the thermodynamic limit. 



the number of grains dissipated at the m-th (last) generation of the n-th avalanche, in 
order to have stationarity, the following balance equation: 

Z{n + 1) = Z{n) + l-<;{p,n), (6) 

or 

p(n + l)=p(n) + i^^^, (7) 

must hold. The average number of grains dissipated at the boundary is (<^(p, n)) = (2p)™ 
(2"^ sites at the boundary, each one occupied with probability p [79]). For each avalanche, 
^{p,n) = (2p)™ -|- ri{p,n), where ri{p,n) is a Gaussian white noise. Plugging this into 
Eq.(I7]) and taking the continuum hmit for n, one can formally write: 

dp _1- (2p)™ r]{p, n) 

whose deterministic part has a stable fixed point at p* = 1/2 = Pc- Accordingly, in 
the thermodynamic limit (in which the effect of fluctuations can be neglected [88]), the 
dynamics attracts p to its critical value (see Fig. |3] left) and the width of the fluctuations 
of p around p = p* = Pc decreases with increasing system size (see Fig. [3] right). This 
simple (conserving) branching process self-organizes to its critical point. 

4-2. Dissipative case 

Borrowing still from Lauritsen, Zapperi and Stanley [90], let us introduce a non- 
vanishing bulk-dissipation rate into the Manna model and, as a consequence, into its 



Self- organization without conservation 



true or just apparent scale-invariance? 15 



self-organized branching process representation. Each offspring (not necessarily in the 
last generation) is removed from the system with probability e; the effective branching 
probability becomes q = p{l — e), and the criticality condition is g = 1/2, or 

and (<^(g, n)) = {2q)™- = (2p(l — e))™'. Eq.(I7]) transforms into: 

p{n + l)=p{n) + — . (10) 

where K{q,n) is the total amount of grains dissipated in the bulk. A simple calculation, 
synthesized in Appendix B, leads to the following equation for the evolution of p 

where 



N N[{l-p)+ pe] 



^-P^(^-^')'""-2(2p(l-e)r 



1 - 2p{l - e) 

(12) 

and, as above, the noise amplitude is N-dependent. After some simple algebra and 
omitting the noise term, Eq. ffTT]) can be rewritten as 

I - (1^(1 - (13) 

with X = 2p{l — e). It is straightforward to check that the only stable fixed point of 
Eq. ([T^ is p* = 1/2 [91J. Therefore, as Pc = 1/2(1 — e), the self-organized dynamics leads 
to a sub-critical point p* = 1/2 < p^, the fixed-point branching ratio is less than unity, 
and the process propagates only for a finite number of generations for any non- vanishing 
value of e [90] • It is only in the conserving limit, e = 0, that the self-organized value and 
the critical point coincide; otherwise there is self- organization to a sub-critical point. 



4-3. Dissipation and loading 

In a recent series of papers [35], Juanico and collaborators introduced a background 
dynamics into the self-organized branching model with dissipation. These authors 
consider a dissipative version of the Manna sandpile rules in the following way: with 
probability a, an active site transfers 2 grains to 2 different randomly chosen neighbors; 
with probability /3, only one grain is transferred to one neighbor while the other one is 
dissipated; finally, with probability e = 1 — a — (3, the two toppling grains are dissipated. 
For simplicity and to easy comparison with the calculations above, we fix /3 = 0; the 
critical branching probability becomes pc = 1/2(1 — e) = l/2a. 

A background dynamics is implemented in ^35j by introducing a rate A for a stable 
site to be turned into a critical one {z = ^ z = 1), and a rate for the opposite 
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Figure 4. Fixed point of the evolution of Eq. p^ in the self-organized branching 
process as a function of A, for a fixed dissipation parameter e = 0.25. Note that only 
for A = Ac(e) « 2.288 10"'', the fixed point corresponds to the critical value of p, 
Pc(e) — 1/2(1 — e) = 2/3. Otherwise, the system self-organizes to a subcritical or to a 
supercritical point. 



transformation. From now on, and without loss of generality, we restrict ourselves to 
the "loading" process (which increases the energy) and fix = 0. Neglecting the noise 
term, it is straightforward to arrive at the following evolution equation for p, 



„m 



with X = 2pa. It has a stable fixed point at 



X = J^l^. (15) 

{l-x)N 1-p ^ ^ 

Note that, at Pc x is equal to 1 and thus, if A is fine-tuned to 

then the fixed point of Eq.(IT4l) becomes pc, i.e. the pair (Ac, Pc), with Ac given by Eq.(IT6ll 
andpc = l/2a, fulfills Eq.f lTSl) and the self-organized branching process becomes critical. 
On the other hand, fixing A > Ac (resp. A < Ac) the fixed point becomes supercritical 
(resp. subcritical) as illustrated in Fig. HI In conclusion, by carefully tuning the loading 
parameter to exactly compensate the effect of dissipation, the system self-organizes to 
its critical point. This process, requiring an explicit parameter tuning, cannot be called 
bona fide self- organization. Finally, note that both, the critical loading parameter Ac and 
the driving rate 1/A^, vanish in the large system size limit (m — )■ 00, — )• c)o), while the 
ratio "loading over driving" (Ac = 771^1-^) diverges. These conditions are analogous to 
those usually imposed to forest-fire models {p, / — 0, p// — )■ 00). The calculation above 
illustrates that such conditions are necessary but not sufficient to achieve criticality in 
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the absence of a conservation law in any finite system: a size- dependent fine tuning is 
also required. 

Even if Juanico and coauthors claim to have designed critical self-organized models, 
all branching processes studied by them [35] are similar in spirit to the example above: 
close inspection of their rules reveals an underlying parameter fine tuning in all the 
different variations they study. For example, for the full model described above (i.e. 
with /3 7^ 7^ 0), Juanico et al. explicitly make the "convenient" choice of parameters 



analogous to Eq.( fT6l) . Actually, in [35], Eq. (JT711 is obtained by fine-tuning the A^- 
independent terms to the critical point, while the last term in the flow equation, 
proportional to is neglected. Not surprisingly, the resulting model converges to its 
critical point in the thermodynamical limit; but, requiring careful tuning, it cannot be 
properly called critical self-organization. 

Summing up, in this section we have illustrated that the conserving self- 
organized branching process shows asymptotically critical dynamics, while its dissipative 
counterpart self-organizes to a subcritical point. Introducing a loading mechanism, 
dissipation can be compensated and criticality restored if and only if the loading rate 
is fine-tuned to a precise size-dependent value. Otherwise, the system self-organizes 
generically either to a subcritical point or to supercritical one. 

5. A full description: Langevin theory of SOC 

Having already studied two different mean-field like approaches, in this section we 
discuss the complete Langevin theory of self-organizing systems. This theory explains 
the origin of the underlying critical point beyond mean-field, its universality in any 
dimension, as well as the key mechanism producing SOC. 

First (Subsection 15. ip . we review the existing absorbing-phase-transition Langevin 
picture of conserving SOC models. Then (Subsect. 15. 2p . we introduce bulk-dissipation 
and extend the theory to non- conserving systems. Finally (Subsect. 15. 2p . a loading 
mechanism is introduced to elucidate the behavior of non- conserving self-organized 
systems. We emphasize the substantial differences with respect to the conserving case. 

5.1. Langevin theory with conservation 

The main idea to construct a stochastic theory of conserving SOC is to "regularize" 
sandpiles (and related systems) by switching off both boundary dissipation and slow 
driving [22l |25l HOI Hi] |92]. In this way, the total amount of sand or "energy", E, in the 
pile becomes a conserved quantity, and can be retained as a control parameter. Indeed, 
in the so-defined fixed- energy ensemble" and for large values of E, the system is in an 
active phase with never-ending relaxation events. Instead, for small values of E it gets 
trapped with certainty into some absorbing state [93] where all dynamics ceases (i.e. all 
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sites are below threshold). Separating these two regimes there is a critical energy, Ec, 
at which an absorbing phase transition takes place. In this way self-organized criticality 
is related to a standard phase transition [22], [911 195] . 

It has been shown [221 SD ES] that such a critical value, coincides with the 
stationary energy density to which the original self-organizing sandpile converges. In 
other words, the energy around which the standard sandpile (i.e. including slow-driving 
and boundary dissipation) fluctuates is the critical point of the "fixed-energy sandpile" . 
Furthermore, the width of fluctuations becomes smaller and smaller with increasing 
system size (see left part of Fig. [5]), guaranteeing that in the thermodynamic limit the 
original sandpile self-organizes to criticality. 

This connection between "driven/dissipative" systems and their "fixed-energy" 
counterparts permits us to relate avalanche exponents to standard critical exponents (see 
[97] for scahng relations) and to rationalize the critical properties of SOC systems from 
the broader point of view of standard non- equilibrium (absorbing-state) phase transition 

[2211111 EH]. 

Using this approach, it has been established that stochastic sandpiles do not 
belong to the robust directed percolation (DP) class, prominent among absorbing phase 
transitions, but to the so-called "conserving-DP" ( C-DP hereafter) or Manna class. This 
class is characterized by the coupling of activity to a static conserved field representing 
the conservation of sandgrains [221 HU |99l 1100] . The field theory or set of mesoscopic 
Langevin equations proposed under phenomenological grounds to describe this class is: 

dtp{x, t) = ap- hp^ + upE + DV'^p + ay/pr]{x, t) 
dtEix,t) = DEV'pix,t) ^ ' 

where p{x,t) is the activity field (characterizing the density of grains above threshold), 
E{x,t) is the locally-conserved energy field, a,b,co, D,a, and De are parameters and 
f]{x,t) is a Gaussian white noise. Some dependences on {x,t) have been omitted to 
unburden the notation. 

Note that in the C-DP class, two fields are required for a Langevin representation: 
the activity field representing grains/energy /force above threshold and the background 
or energy field describing the local amount of grains/energy /force. Nevertheless, it is 
possible to stick to a single field description by integrating out the energy equation. 
This generates two extra terms for the activity equation 

u}pE{x,0)-up [ dt'V'^pit'). (19) 
Jo 

The second, history dependent (non-Markovian) term describes the tendency of sites 
that have been less active than their neighbors in the past to be more susceptible for 
activation (e.g. in a sandpile, if the neighbors of a given site have toppled, the site is 
very likely to overcome the threshold). The single equation for the activity reads 

dtp{x, t) = (a + uE{x, 0))p - bp^ - up dtV^p{x, t') + DV^p + (r^r]{x, t); 

(20) 
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non-Markovianity is the price to pay for removing the energy field |1U1] . 

The C-DP class described either by Eq. ffTSj) or by Eg. (1201) has a critical dimension 
dc = 4: and embraces not only stochastic sandpiles, but also (among other examples) 
some conserving reaction-diffusion systems, for which the equations above can be 
explicitly derived from the microscopic dynamics [991 HOP] ■ 

Eq. ffTSl) can be studied either (z) in spatially extended systems, (ii) using random 
neighbors, or (Hi) in a globally, all-to-all, coupled version in which the Laplacian is 
replaced by p — p{x, t) ~ ~ Pi^y ^)]- These last two are useful to construct 

mean-field approximations of the full (spatially extended) theory. 

5.1.1. Relation with other universality classes First, note that the well-known theory 
for directed percolation (i.e. the Reggeon field theory [93]) is recovered upon fixing 
w = in Eq. ffTSj) . The additional conserved field turns out to be a relevant perturbation 
altering the critical behavior of systems in the DP class [221 SD HOOj - 

On the other hand, for the sake of completeness, let us just briefiy mention that 
the C-DP class is fully equivalent to the pinning/depinning transition of interfaces in 
random media, i.e. the Quenched Edwards Wilkinson [231 11021 I103j . The absorbing 
(resp. active) phase maps into the pinned (resp. depinned) one. Exploiting the mapping 
between these two descriptions (of a unique underlying physics), critical exponents for 
Eq. (fT8!) can be deduced from existing renormalization group results for interfaces [104j . 

It is worth stressing that, in terms of interfacial models, conservation of energy is 
equivalent to interface translation invariance; for instance, a dissipative term like —ep 
introduced in Eq. ffTSl) would map into a term —eh (where h = h{x,t) is the interface 
height) in the Quenched Edwards Wilkinson equation, which breaks such an invariance. 

In this respect, recent experimental evidence of self-organized critical behavior 
(including finite size scaling), obtained for avalanches in type II superconductors [20] , 
give critical exponents compatible with those of the C-DP class. Barkhausen noise [105] 
and acoustic emission in fracture [106j are other related examples. 

5.1.2. Conserved SOC as an absorbing state transition. Within this framework, the 
way self-organized criticality works is as follows (see left part of Fig. ^ [221 E] • if 
sandpile is in its absorbing phase {E < E^., where E stands for the spatially averaged 
value of E{x) in the steady state) then, owing to the driving mechanism, the energy 
is slowly increased until, eventually, the active phase {E > E^) is reached. At this 
point avalanches are triggered and they re-structure the sandpile energy configuration. 
Avalanches may dissipate energy at the open boundaries, until eventually the system 
falls back into an absorbing state, the avalanche stops, and slow driving acts again 
restarting the cycle. In this way, the sandpile is expected to fiuctuate around its critical 
point, Ec, with excursions to either the active or the absorbing phase as sketched in the 
upper-left part of Fig. [5l 

To have a numerical confirmation of this, Eq. (fT8l) can be interpreted as in SOC, i.e. 
one can implement slow driving and dissipation at infinitely separated time-scales, and 
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integrate the equation using the efficient algorithm introduced in |107j . In particular, 
one considers an absorbing configuration (p(x, t = 0) = 0) and open boundaries, then 
add a small amount of activity /energy, 6, to a given site, xq'- 

p{xo, t = 0) ^6, E{xo, t = 0) ^ E{xo, t = 0) + 6; (21) 

this generates an avalanche, which evolves according to Eq.( |T8i) . Iterating this process, 
one obtains the distribution of values E sampled during avalanches, shown in the lower- 
left panel of Fig. [5l Observe that, as dissipation and driving become arbitrarily small 
by increasing system size (actually, they are infinitesimally small in the thermodynamic 
limit), the degree of penetration into the active and absorbing phases is arbitrarily small, 
the distribution of E becomes more and more peaked, and the system is arbitrarily 
close to its critical point. Moreover, the avalanche exponents measured by means of 
numerical integration of Eq. fllSI) at such a steady state coincide with (or can be related 
to) those obtained by performing standard fixed-energy simulations of Eq. fllSI) at its 
critical point [1071 1103] . Note the obvious analogies between this picture and the self- 
organized branching process described above: E is the equivalent of p, i.e. the self- 
organized control parameter; the critical point Ec corresponds to the critical branching 
probability pc- 

The advantage of Eq. ffTSl) as a theory for SOC with respect to the self-organized 
branching process is that, while this last is a mean-field theory explaining qualitatively 
self-organization but failing to justify critical exponents in spatially extended systems, 
Eq. lfTSl) is a full theory including fiuctuations and spatial-dimensionality. It provides 
accurate estimates for avalanche exponents in any dimension and opens the door to 
field theoretical analyses. Furthermore, Eq. fllSI) . considered on a random- neighbor or 
an all-to-all coupling, constitutes a sound mean-field description of conserving SOC, 
equivalent to those in the preceeding sections. 

To end up, note, once again, the essential role played by conservation in this theory. 
The underlying phase diagram sketched in the left part of Fig. [5] relies on the averaged 
energy being a control parameter. If there was a non-vanishing bulk dissipation, the 
energy would change continuously during avalanche evolution. In the next sub-section, 
we shall explore how this affects the absorbing-phase-transition picture of SOC. 



5.2. Langevin theory with bulk- dissipation 

To tackle the problem of non-conservation within the absorbing-state Langevin 
framework, we need to modify Eq.f JTSl) to allow for bulk dissipation. Introducing in 
Eq. ffTSl) the leading dissipative term, —ep{x,t), and neglecting higher order corrections, 
the resulting set of equations becomes 

dtp{x, t) = ap- hp^ + ujpE{x, t) + DV'^p + ay/pr]{x, t) 
dtE{x,t) = DEV^pix,t)-ep{x,t) ^ ' 

which is, obviously, non- conserving owing to the activity-dependent energy leakage. 
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Figure 5. (Upper part) Sketched diagrams illustrating the mechanism of self- 
organization for conserving dynamics (left panels) and for non-conserving dynamics 
with loading (right panels). Left-up: In the "fixed-energy" ensemble, there is an 
active and an absorbing phase, separated by a critical point, Ec- Slow driving and 
dissipation make sandpiles (and related conserving systems) fluctuate around their 
associated fixed-energy counterpart critical point. Left-down: the distribution of E 
during avalanches is plot for simulations of the conserving Langevin theory, Eq. ()18|) . 
using various system sizes (2^^, 2^^, and 2^"^). Other parameter values are a = 0.423, 
b = uj = l, D = D2= 0.25, a — \/2; the distributions become progressively peaked 
around E^ upon increasing system size (the value of E has been normalized with 
Ec for each size). Right: Plots analogous to their counterparts to the left, but for 
non-conserving systems. Right-up: In this case, there is no active phase, but just a 
phase in which perturbations/avalanches spread. Results in right-down correspond 
to simulations of Eq. ([22|) . using Eg. ([25]) and Eq. ([26|) . with Emax = 1-3 (rest of 
parameters as above) . Observe the much broader distributions appearing in this case: 
the system hovers around the spreading critical point with large excursions into both, 
the supercritical and the subcritical phase. 
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Integrating in time the equation for the energy field, the following extra terms for 
the activity equation are generated: 

up[E{x,0) + [ dt'V'^p{x,t') -e [ dt' p{x,t')]. (23) 
Jo Jo 

The second term, dominant in Eq. fl20|) . becomes a higher order correction here, i.e. 
it is irrelevant in the renormalization group sense as compared with the third, non- 
Markovian, term. The last one is well known to be the leading non-linearity in the 
dynamical percolation universality class [108] . From this perspective, it is no wonder 
that the critical behavior of some non- conserving SOC models (e.g. forest-fires) has 
been related to (dynamical) percolation in the literature jl09]. Such a class, whose full 
(one-field) Langevin equation is: 

dtp{x, t) = {a + ujE{x, 0))p - hp^ - eup [ dt'p{t') + DV^p + a^T]{x, t), (24) 

J 

describes the spreading properties of epidemics with immunization, etching of disordered 
solids |110] . and some aspects of spreading in systems with many absorbing states 
without a conservation law (see for more details). The term —ecop j^dt' p{t') ensues 
that regions already visited by activity become less prompt to be active in the future. 
Owing to this term, Eq.( l22l) cannot sustain an active phase. However, even if it lacks a 
stable active phase, it exhibits a spreading phase transition separating a phase in which 
seeds of activity propagate indefinitely (in the form of rings of expanding activity; i.e. 
defining an "annular growth" phase) from an absorbing phase in which they do not [108] 
(see the diagram at the right part of Fig. [5l). 

This spreading transition, whose critical dimension is d^ = 6, is controlled by 
the initial state in which the seed of activity is placed; indeed, the initial energy, 
uE{x,t = 0)p{x,t), is a mass term in Eq.f l23p . If the initial value of E is large enough 
(i.e. favorable environment), then avalanches tend to propagate, while if it is small, they 
do not. Separating these two regimes there is a critical point for spreading propagation 
at some value E*{x,0). 

In summary, the introduction of a non-vanishing dissipation rate affects in a 
relevant way the critical behavior of self- organizing systems; depending on the initial 
condition, the dissipative Eq. fl22|) can be in the propagating/supercritical or in the non- 
propagating/subcritical phase of a dynamical percolation phase transition. 

5.3. Full theory: dissipation and loading 

Now, we are in a good position to understand in depth the role of the loading mechanism 
in dissipative models of SOC, within the Langevin framework. To do so, let us 
complement Eq. (!22l) with a specific prescription on how to change the background energy 
field between avalanches, i.e. let us consider a loading rule, as for instance: 



E{X, t = 0)^ E{X, t = 0)+ -f{Ema. - E) 



(25) 
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Figure 6. Histograms of E, sampled during avalanches, in a numerical integration 
of Eq. (f22|) for (left) e = 0.1 and different values of the loading parameter 7 (as soon 
as 7 > 7c the distribution develops a double-peak shape) and for (right) 7 = 0.1 and 
different values of e. Other parameter values are as in Fig. [5l 

(where E^ax is a parameter and E the average energy in the system) and a driving rule: 

p{xo,t = 0)^6. (26) 

Eq.( l25|) and Eq.f l26p define one possible loading mechanism for Eq. fl22|) : other choices 
are, of course, possible (for instance, the loading mechanism could also act "during" 
avalanches). The results presented in what follows are generic, essentially independent 
of such a choice. 

As shown above, Eq.( l22l) can sustain avalanches propagating indefinitely (up to 
system size), provided that the initial energy is large enough. Therefore, considering 
large values of 7 or of Emax in Eq.( l25|) ). the system becomes supercritical for avalanche 
propagation. Instead, small initial densities lead to subcritical propagation. 

To illustrate this and the forthcoming discussion, we have performed computer 
simulations of Eq.(l22l) using the parameter values specified in the caption of Fig. O and 
complemented with the loading and driving rules Eq.f l25p and Eq. fl26|) . For the sake 
of simplicity, we have considered sites with an all-to-all (mean-field) coupling. To 
check the robustness of our conclusions we have also studied a random-neighbor version, 
obtaining very similar results (not shown). As before, the equation has been integrated 
using the algorithm introduced in jl07]. 

For a dissipation parameter e = 0.1 and fixing E^ax = 1-3, we find a critical point 
for spreading at some value of 7, 7c ~ 2 x 10^^, which generates, on average, an initial 
energy E* ^ 1.277(5). At such a critical value, power-laws for avalanche and spreading 
exponents are obtained. In complete analogy with Fig. [H for smaller values of 7 (and 
hence, smaller values of the initial average energy density) the avalanche size distribution 
has an exponential cut-off (sub-critical), while for larger values the distribution develops 
a bump for large sizes (super-critical) (results not shown). 

We have constructed histograms of the average energy by sampling E during 
avalanches in computer simulations (see Fig. [6]). Typically, for short times the system 
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is in the right side of the distribution and, as the avalanche proceeds and dissipation 
acts, E moves progressively leftward. This shifting generates a broad distribution of 
energy values for any system size. After the avalanche stops, the system is "loaded" 
again (Eg. (125!) ). a new avalanche is triggered (Eq.( l26|) ). and so on. In this way, the 
system is kept hovering around the critical point, with a broad distribution of values, 
as illustrated in Fig. [5] and Fig. [61 

In Fig. [6] (left panel), histograms for e = 0.1 and various values of 7 are plot. As 
long as the loading is sufficiently strong (7 > 7c ~ 2 x 10^^) the distribution develops a 
double-peak structure overlapping with both the propagating and the subcritical phases. 
Instead, for smaller values of 7 (7 = 10~^, in the figure), the loading is too small, and the 
histogram overlaps only with the subcritical phase. In Fig. [6], (right panel) histograms 
for 7 = 0.1 and various values of e are plot; in all cases, there is a double-peak structure. 
Observe that, for large dissipation rates, the system gets deeper into the absorbing phase. 

It is important to remark, that (as illustrated in Fig. [5]) the dynamics is rather 
different from its conserving counterpart: while, in the conserving case, fluctuations 
around Ec decrease in amplitude with system-size (Fig. [5l left-panels), in the non- 
conserving case the histograms remain broad, even in the thermodynamical limit (see 
Fig. [5l right panels); i.e. large variations around the critical spreading point persist for 
any system size. 

We caution the reader that, within an avalanche, the process is not stationary 
(energy decreases) and, therefore, the histograms shown in Fig. [5] and Fig. [6] cannot be 
properly interpreted as probability distribution functions. 

For this reason, we have also constructed stationary (steady state) histograms 
for (i) the distribution of the average initial energy for avalanches and (ii) the 
distribution of E values after avalanches. Any of these can be used, as well, to 
illustrate the differences with the conserving case. For example, Fig|7] shows that the 
background in which avalanches are started is generically non-critical, but is broadly 
distributed around the critical point. Using this information, one can make the educated 
guess that the associated avalanche size distribution (or any other quantity measured 
for avalanches/spreading), with such a distribution of initial conditions, will be a 
convolution of different subcritical and supercritical curves weighted with the above 
distribution of initial energies. Let us emphasize the lack of any mechanism tuning the 
system to criticality: the energy is initially set to some arbitrary value (controlled by 
the parameter 7). If and only if the initial density is fine tuned to the critical value for 
spreading of activity, E*{x, 0), the system is at the critical point for avalanche spreading. 
Otherwise, for larger values of 7 the system is initially supercritical, while for smaller 
values it is sub-critical. 

In conclusion, flne-tuning of the loading mechanism is required to have critical 
spreading in non- conserving systems, and it is controlled by a dynamical percolation 
critical point. The important point to stress is that, even if the initial condition 
is not critical for generic values of 7 and Emax, the "hovering- around-the- critical- 
point" mechanism, illustrated in Fig. [5] (right panels), keeps the dynamics effectively 
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Figure 7. Histograms showing the distribution of initial conditions, E, for avalanche 
propagation for various system sizes in numerical simulations of Eq. (|22p endowed with 
the loading mechanism Ea. ((25|) (parameter values are as in FiglS]). Observe that the 
distribution remains broad in the thermodynamic limit. 

not far from criticality, but not at criticality, for a broad range of parameter values. 
Furthermore, contrarily to the conserving case (Fig. [5l left panels), the energy 
histograms do not tend to a delta-peak function for large system sizes. Large excursions 
into both the supercritical and the subcritical phases persist in the thermodynamic 
limit. 

5.3.1. Revisiting the mean-field Pruessner- Jensen model. Using this insight, we can 
now tackle the open question: how does l/6c{N) scales at criticality in the model 
discussed in Section ISp 

As l/6c{N) is a background energy, it corresponds to a mass term in the Langevin 
equation for the activity, Eq.(l23l). Therefore, to preserve scale-invariance, it needs to be 
scaled with iV as a distance to the critical point and, therefore: 

l/e,iN) ~ (27) 

where ^ is the correlation length and u the correlation length exponent, whose mean- field 
value in the dynamical percolation class is z/ = 1/2 [97] . 

At the upper critical dimension, dc = Q (where hyper-scaling mean-field relations 
are expected to hold), ^ is limited by the system linear size, N^^'^ = N^^^, where N is 
the volume (i.e. total number of sites): 

^ ~ N^/^. (28) 
Putting together these two last equations, we obtain the mean- field scaling result. 




(29) 
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This is in excellent agreement with the empirical result reported in Eq. . Analogously, 
the cut-off of the avalanche size distribution, Sc, needs to scale as 

s, ~ ~ N^f^'^^ = N^/\ (30) 

where we have used the fractal dimension, Df = 4, for mean-field dynamical percolation 
[97] . Again, this prediction is in perfect agreement with the numerical results Eq.(l5]). 

In conclusion, the loading parameter needs to be fine-tuned for each finite size to 
have true scale-invariance; its scaling is inherited from a dynamical percolation critical 
point. 

5.3.2. Spatially extended systems. The numerical results reported in this section 
correspond, as already stressed, to an all-to-all (as well as random-neighbors) coupling 
in Eq. (!22l) . On the other hand, qualitatively similar results can be obtained for spatially 
extended systems: i.e. a broad distribution of the spatially-averaged energy, hovering 
around the critical point. The main difference is that, obviously, the background energy 
becomes heterogeneously distributed in space and, therefore, the situation becomes 
much more involved. 

Eq. fl22l) generates spontaneously regions with higher and with lower values of 
E, which have different propensities to activity propagation: there are locally super- 
critical and sub-critical regions. Patches where avalanches have passed are typically 
less likely to propagate new activity owing to the non-Markovian term in Eq.( l22l) . 
The size distribution of such patches is, accordingly, inherited from the avalanche size 
distribution, creating a complex landscape for further avalanche propagation. This 
scenario is, of course, more complex than in the mean-field one discussed above, but 
the essence of the described phenomenology remains un-altered: the (local) control 
parameter hovers- around a dynamical percolation critical point, with fluctuations that 
do not vanish in the large system-size limit. 

Observe that, in order to tune the system to criticality (as done in the mean- field 
case) one would need in this case to define a more complicated loading mechanism which 
should get rid of the dynamically generated heterogeneities, leading to a homogeneous 
initial-energy state, tuned exactly to its critical value. Using the language of [22], 
one needs to "hire a babysitter" (or a "gardener" using the forest-fire terminology 
[95] ) to keep the spatially extended system sitting (everywhere) at criticality. Once 
such an efficient babysitter is at work, the initial condition is always at the (dynamical 
percolation) critical point in any dimension. 

From this perspective, our theory provides additional support for the claim in [86] 
that the critical density of trees in a forest fire model should coincide with the percolation 
critical density, i.e. in order to observe critical propagation in the forest-fire model one 
should tune the initial background (number of trees) to the corresponding dynamical 
percolation critical density. 

As already pointed out by Grassberger some time ago [86j, in the absence of an 
efficient gardener taking care of local fine-tuning, partial power-laws are still observed 
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in the forest-fire model. This is due to the existence of patches with different densities 
of trees, which appear with a broad spectrum of sizes. Each of such patches hes at a 
different distance of the critical point. The convolution of avalanches propagating in 
such a variety of initial conditions originates a complex pseudo-scaling picture which, 
obviously, does not correspond to strict criticality. A similar picture applies also to the 
(open-boundaries) OFC earthquake model [72], and to self-organized models of neural 
activity (as will be illustrated in a forthcoming paper). 

Summing up: In this section, we have reviewed the standard absorbing state phase 
transition picture of SOC, underlining the special role played by conservation. Then, we 
have introduced a bulk dissipation term and illustrated that the active phase disappears 
and the universality class of the involved spreading phase transition is changed from 
C-DP to dynamical percolation. Dissipation needs to be compensated by a loading 
mechanism (which controls the initial conditions in which avalanches are started) to 
keep the energy balance. If and only if the loading mechanism is perfectly fine-tuned 
to generate a precise initial energy density, true criticality is observed. Otherwise, the 
system just hovers around a dynamical percolation critical point, with large excursions 
into the propagating and the absorbing phases. Contrarily to the conserving case, 
such fluctuations do not disappear in the thermodynamic limit. Strictly speaking, this 
mechanism of self-organization cannot be called critical, we propose to refer to it as 
self- organized quasi- criticality (SOqC). Last but not least, our approach provides a way 
to rationalize the finite size scaling properties of non-conserving self-organized systems, 
as earthquake or forest-fire models, and has allowed us to derive a number of previously 
unknown scaling relations. 

6. Concluding Remarks 

We have shown, by using different levels of description, that non-conserving models of 
self-organized criticality, as earthquake and fire-forest models, are not truly critical. 

First, we have studied a simple mean-field theory, based on an energy balance 
equation, for a non- conserving model introduced by Pruessner and Jensen. It permits 
us to illustrate that, even if one can construct non- conserving models that seem critical 
in the thermodynamic limit, there is no systematic way to have a coherent finite-size 
scaling description of them: a precise fine tuning is required for any finite size to observe 
criticality and to approach the thermodynamic limit in a scale-invariant way. 

Second, we have revisited the mapping of high-dimensional (mean- field) avalanching 
systems into a self-organized branching process, introduced by Zapperi et al. some years 
ago. The underlying idea is that in high dimensions avalanches do not visit twice a given 
site and they can be described as a branching process, whose branching probability 
depends on the energy background. This allows us to write an evolution equation 
for the branching probability for models with slow driving and dissipation. While in 
the conserving case the branching probability converges to its critical value, in the 
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presence of bulk-dissipation the convergence is towards a subcritical point. Introducing 
a loading mechanism (which mimics the growth of new trees in forest fire automata 
or the continuous building-up of stress in earthquake models), we have shown that the 
fixed-point towards which the system self-organizes can be either critical, subcritical, or 
supercritical. Contrarily to previous claims, a fine-tuning of the loading mechanism is 
required to reach criticality within this approach. 

Third, we have introduced a full stochastic description of SOC systems in terms 
of Langevin equations. We have reviewed how conserving systems self-organize to a 
critical point with well known critical exponents (in the conserved-directed-percolation 
universality class). Instead, for dissipative systems with a loading term, the dynamics is 
found to hover around a critical point, with large excursions into the absorbing and the 
propagating phases, which do not disappear in the large system size limit. Therefore, 
such systems are not generically critical. Still, some traits of the underlying dynamical 
percolation critical point can be observed, depending on the loading parameter and 
system-size, i.e depending on how large the excursions into the absorbing/propagating 
phase are. 

All these three approaches provide overwhelming evidence that conserving- 
dynamics is a necessary condition to observe self-organization to criticality. Instead, 
in non-conserving (dissipative) systems equipped with a loading mechanism, a fine 
tuning of a loading parameter is required to have the system sitting at a critical point. 
Otherwise, the system just hovers around a critical point, with broadly distributed 
fluctuations which do not dissappear in the thermodynamic limit: for a broad range of 
parameters, non-conserving systems can be fluctuating in the vicinity of a critical point, 
but not at the critical point. We propose to call this: self- organized quasi- criticality.. 

This conclusion extends to some other non-conserving models of SOC as those for 
synchronization of integrate-and-fire oscillators described in [611 1112[ 1113] , the model 
in |114] . or the model of neural avalanches in [36] (as will be explained in a separate 
publication). 

Is it sensible to refer to self- organized dissipative systems as "critical"? 

The answer to this question is mostly a matter of taste, and depends on what one 
wants to define as criticality. Being strict and calling critical only to systems sitting 
at a critical point (allowing at most for fluctuations that vanish in the thermodynamic 
limit), then non- conserving systems are not truly critical models. 

Being more permissive, one could accept, in principle, the term "critical" to refer to 
systems hovering around a critical point (with persistent excursions into the subcritical 
and supercritical regimes), which exhibit "dirty scaling". 

In order to avoid miss-understandings and miss-conceptions, we strongly favor the 
use of an alternative terminology like "almost criticality" [67], "pseudo-criticality" , or, 
as we said, self-organized quasi-criticality (SOqC) to refer to non-conserving self- 
organized systems, and suggest to restrict the term "critical" for truly scale-invariant 



Self- organization without conservation : true or just apparent scale-invariance? 29 
systems. 

Actually, in many cases, strict criticality might not be required to explain 
empirical (truncated) power-laws distributions observed in the real world and, therefore, 
"self-organized quasi-criticality" remains a useful concept, despite of the somehow 
inappropriate and certainly confusing use (and abuse) of the word "critical" in the 
literature. Under this light, one could reconsider the empirical observations discussed 
in the introduction, as well as similar ones for which power law distributions have been 
reported. A critical inspection of them reveals in many cases that empirical data are 
better described by truncated power-laws rather than by pure power-laws [1]. 



Appendix A: Basic Models in a nutshell 

For the sake of completeness, in this Appendix we present some of the basic toy 
models of self-organized criticality. All of them are defined in a 2-dimensional 
lattice (generalizations to li- dimensions, to random-neighbor or all-to-all couplings are 
straightforward) . 

Sandpile-like Models 

Consider a (height or "energy") variable z{i,j), which takes integer (non-negative) values 
at each site (i,j) of a two-dimensional lattice. The ingredients of sandpiles, ricepiles, 
and related models can be sketched as: 

• Slow Driving: A small input of energy is externally introduced into the system 
(grains are dropped), usually at a single site: z{i,j) — )■ z{i,j) + 1. 

• Activation: A site receiving energy, stores it until a given threshold, zthr, is 
exceeded and the site is declared active; otherwise nothing happens and the system 
is driven again. 

• Relaxation (or "toppling"): Each active site redistributes (all or a fraction of) its 
accumulated energy among its neighbors. A relaxation event can trigger a chain 
reaction or avalanche by activating its neighbors and so forth. 

• Boundary Dissipation: When redistribution events reach the (open) boundaries of 
the system, energy is dropped off. 

• Iteration: When activity has ceased, the avalanche stops, and a new external input 
is added. The driving/dissipation cycle is iterated until a statistically stationary 
state is reached. 

For each specific model, the relaxation rules and some other details can change, giving 
rise to a zoo of models. We enumerate some of the more commonly studied ones: 

The Deterministic Bak-Tang-Wiesenfeld Sandpile Model: The threshold is 
fixed to Zc = 4: {2d in d-dimensions) . Active sites relax according to: 
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z{i',f)^ z{i',f) + l 



for all (i', j') n.n. of («, j), 



(31) 



that is, an active site is emptied and its grains are deterministically redistributed 
amongst its nearest neighbors. The energy at sites out of the system is fixed to 2; = 0, 
enforcing boundary dissipation when boundary sites topple. Avalanches measured in 
such a steady state were originally claimed to be critical [6]. Later work showed that 
actually, owing to the deterministic nature of the model and, as a consequence, to the 
existence of many toppling invariants and breakdown of ergodicity [115j the system is 
not truly critical but exhibits anomalous multi-scaling |116j . Other authors (see, for 
instance, |117j ) suggested that avalanches do not obey any type of scaling whatsoever. 
To avoid the pathologies associated with deterministic rules, we focus all along this 
paper on stochastic models. 

The Stochastic Manna Model: The dynamics of the Manna model is similar 
to the deterministic BTW, but the threshold is fixed to = 2 in any dimension and 
stochasticity is introduced in the redistribution rule [H]: 



This is the "inclusive" version of the model. In its "exclusive" version, the two grains 
are forced to go to different neighbors. Both of these versions self-organize the model 
to a critical point in the C-DP (or "Manna" ) class in any dimension p3l [99| 1100] . 

The Stochastic Oslo Model: This (ricepile) model has annealed random 
thresholds at each site: every time a site topples, a new threshold Zthrihj) = 2, 3 is 
randomly chosen with equal probabilities [15] . Redistribution of energy units (grains) is 
done as in the Manna model, in a stochastic fashion. These rules lead to self-organization 
to a critical point in the C-DP universality class with rather clean scaling [151 EH HOOj . 

The Olami-Feder-Christensen Earthquake Model 

A continuous "force" or "energy" F{i,j), randomly distributed between and a 
threshold value, Fthr, is initially assigned to each site {i,j) of a two-dimensional lattice. 
During the driving step, the site with maximum force, F^"-^, is identified, and the force 
of all sites is increased by Fthr — ^Ij""^, creating (at least) one seed of activity (this 
is equivalent to advancing all sites at a fixed constant velocity until one of them, the 
maximum, reaches the threshold). Contrarily to sandpile models, this driving affects all 
sites. An active site relaxes according to: 



where a < 1/4 is a bulk- dissipation parameter. The model is conserving only for a = 1/4. 
The force at sites out of the system is fixed to F = 0, entailing boundary dissipation. 



z{i,j)-2 
z{i',f)^ z{t',f) + l 



for 2 randomly chosen {i',j') n.n. of {i,j). 



(32) 






F{t',f) + aF{t,j) V {t',f) n.n. of (2, j), 



(33) 
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The stationary state of such a system was claimed to be critical for a broad range of 
values of a [H] , but recent analyses disprove such a claim (see for instance [551 112] and 
Section [1]). 

The Drossel-Schwabl Forest Fire Model 

Sites of a two-dimensional lattice can be either empty (z = 0), occupied by a tree 
{z = 1), or by a tree on fire [z = 2) [50]. The dynamics proceeds as follows: 

• At every empty site, a tree grows with probability p, and the site becomes 
"occupied" : z = z = 1 

• Initial spark: A tree not surrounded by any fire becomes a "burning" tree (e.g. 
lightening) with probability /: z = 1 ^ z = 2 

• A burning tree sets on fire all its nearest neighbors, and it becomes empty: 
z = 2^ z = 

• The "fire avalanche" proceeds burning all trees in contact with fires. 

• When fire ceases the dynamical processes is re-started. 

The relevant parameter is //p which sometimes is called 1/6 [83]. In the implementation 
of the model that we use, 1/^ is the number of trees grown between two consecutive 
ignitions. Despite of the initial claims and various forms of reported "anomalous 
scaling" , the most recent studies revealed absence of generic scale-invariance [HHl [86] . 



Appendix B 

To compute the mean number of grains lost in the bulk in the dissipative self-organized 
branching process [HE], let us first consider the number of off-springs, ip, not-occupied 
owing either to absence of branching or to bulk-dissipation: 



m—1 




(34) 



then, the average fraction corresponding to bulk-dissipation is 




(35) 



Noting that the avalanche size is s = Y2j 



( m-l 



m 




< 



(36) 



m—1 



m 




\ k=0 



1=1 
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then (fi = 1 + s — 2^m. Using the definition of s: 



k 1 - (2p(l - e)) 



m+l 



k=0 



k=0 



1 - 2p(l - e) 



Plugging Eg. (1371) into the equation for (p 

and, from this and Eq. (l35|) : 

{K{q,n)) -- 



i + i^Mizi)r:_2(2p(i-.)r 



1 - 2p(l - e) 



pe 



[l-p) +pe 



(37) 



(38) 



(39) 



leading, once the continuum limit for n has been taken and fiuctuations included, to 
Eq.(IIID. 
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